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Bias and the power spectrum beyond the turn-over 

Ruth Durrer 1,2 , Andrea Gabrielli 3,4,5 , Michael Joyce 6,7 and Francesco Sylos Labini 8,9 

ABSTRACT 

Threshold biasing of a Gaussian random held gives a linear amplification 
of the reduced two point correlation function at large distances. We show 
that for standard cosmological models this does not translate into a linear 
amplification of the power spectrum (PS) neither at small k not at large k. 

For standard CDM type models the “turn-over” at small k of the original 
PS disappears in the PS of the biased held for the physically relevant range 
of threshold parameters v. In real space this difference is manifest in the 
asymptotic behaviour of the normalised mass variance in spheres of radius 
R , which changes from the “super-homogeneous” behaviour cr 1 2 (R) ~ to 
a Poisson-like behaviour cr 2 (R) ~ R~ 3 . This qualitative change results from 
the intrinsic stochasticity of the threshold sampling. While our quantitative 
results are specihc to the simplest threshold biasing model, we argue that our 
qualitative conclusions should be valid generically for any biasing mechanism 
involving a scale-dependent amplification of the correlation function. One 
implication is that the real-space correlation function will be a better instrument 
to probe for the underlying Harrison Zeldovich spectrum in the distribution of 
visible matter, as the characteristic asymptotic negative power-law £(r) ~ — r ^ 4 5 6 7 8 9 
tail is undistorted by biasing. 
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Subject headings: galaxies: general; galaxies: statistics, biasing; cosmology: 
large-scale structure of the universe 

The concept of bias has been introduced by Kaiser (1984), primarily to explain the 
observed difference in amplitude between the correlation function of galaxies and that of 
galaxy clusters. In this context the underlying distribution of dark matter is treated as a 
correlated Gaussian density held. The galaxies of different luminosities or galaxy clusters, 
are interpreted as the peaks of the matter distribution, which have collapsed by gravitational 
clustering. Different kind of objects are selected as peaks above a given threshold, with 
a change in the threshold selecting different regions of the underlying Gaussian held, 
corresponding to huctuations of differing amplitudes. The reduced two-point correlation 
function of the selected objects is then that of the peaks ^(r), which is enhanced with 
respect to that of the underlying density held £(r). In a previous paper (Gabrielli, Sylos 
Labini & Durrer, 2000 - GSLD00) some of us have discussed the problematic aspects of this 
mechanism. In particular the amplification of the correlation function is in fact only linear 
in the regime in which £„(r) <C 1 (see also Politzer & Wise, 1984). In the region of most 
observational relevance (where ^(r) 3> 1) the correlation function is actually distorted at 
least exponentially. Furthermore we have drawn attention to the fact that the amplification 
of the correlation function by biasing rehects simply that the distribution of peaks is more 
clustered because peaks are exponentially sparser. 

In this letter we discuss a different aspect of this model for bias. We are interested 
in understanding the effect of biasing on the power spectrum (PS). In particular we 
address here a qualitative change that is caused to the matter perturbations in in 
standard cosmological models. In real space this change manifests itself in a change from 
sub-Poissonian behaviour of the mass variance at large scales in the underlying density 
held, to Poissonian behaviour of the same quantity for the biased held. In k space this 
implies a distortion of the PS at small k. Our analysis shows that this effect can be very 
important observationally, as it can make the “turn-over” in the dark matter PS disappear 
from the PS of visible matter. Furthermore it shows the importance of measuring not just 
not just the PS of visible objects, but also their real space correlation properties. Earlier 
works about the effect of biasing in the PS can be found e.g. Coles (1993) and Scherrer & 
Weinberg (1998). 

Before considering threshold biasing, we recall the relevant part of the analysis given in 
Gabrielli, Joyce & Sylos Labini (2002) (hereafter GJSL02). In this paper we have discussed 
the meaning of the condition P(0) = 0 satisfied by the PS of all current standard type 
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cosmological models (with Harrison-Zeldovich like spectra P(k) ~ k at small k). While 
this point is often noted in the cosmological literature (see e.g. Padmanabhan 1993), its 
significance and implications are not correctly appreciated (see GJSL02 for discussion). 

It implies the requirement that the integral over all space of the correlation function 
vanishes, meaning that in the system there is an exact balance between correlations and 
anti-correlations at all scales. This is a highly non-trivial, non-local, condition on the 
distribution. Its specificity can be highlighted by the following classification of all stationary 
stochastic processes into three categories: (i) For -P(O) = oo the fluctuations are like 
those in a critical long range correlated system, (ii) for -P(O) = constant > 0 the system 
is Poisson-like at large scales e.g. any short-range positively correlated system such as a 
quasi-ideal gas at thermal equilibrium, and (iii) for -P(O) = 0 the system is what we have 
termed “super-homogeneous”. The reason for the use of this last term comes from the fact 
that the three categories are distinguished most strikingly in real space by the large distance 
behaviour of the mass variance in spheres, as one can show that P(0) = limy^oo (( AJV W)) ) 
where ((AM(h)) 2 ) is the mass variance in a volume V (and p Q the mean mass density). In 
the Poisson type distribution this variance is proportional to the volume of the sphere, while 
in the first category (critical systems) it grows more rapidly (with a limiting behaviour of the 
volume squared), while in the last (super-homogeneous distributions) the growth is slower 
than the Poissonian one. In particular the case of the H-Z spectrum marks the transition 
to the limiting slowest possible growth of this quantity for any stochastic distribution of 
points (Beck 1987), which is a growth proportional to the surface of the sphere. 


These super-homogeneous distributions are encountered in various contexts in 
statistical physics. They are described in this context as glass-like: they are highly ordered 
distributions like a lattice, but with full statistical isotropy and homogeneity. In Gabriclli 
et al. (2002) an example of a system with such correlations at thermal equilibrium is given, 
and a modification of this same system which should give precisely the correlation of a 
standard cosmological model is described. 


We now turn to the threshold biasing mechanism. Following Kaiser (1984) we consider a 
stationary, isotropic and correlated continuous Gaussian random held, <5(x), with zero mean 
and variance a 2 = (A(x) 2 ) in a volume V as V —> oo. The marginal one-point probability 

5 2 

density function of 5 is V(S) = ■ Using V(5), we can calculate the fraction of 

the volume V with <5(x) > ua, Q\{v) = /“ V(5)d5. The correlation function between the 
values of h(x) in two points separated by a distance r is given by £(r) = (h(x)A(x + rn)). 
By definition, £(0) = a 2 . In this context, stationarity means that the variance, cr 2 , and the 
correlation function, £(r), do not depend on x. Statistical isotropy means that £(r) does 
not depend on the direction n. The goal is to determine the correlation function of local 
maxima from the correlation function of the underlying density held. The problem can be 



4 


simplified (Kaiser 1984) by computing the correlation of regions above a certain threshold 
va instead of the correlations of maxima. However, these quantities are closely related for 
values of v significantly larger than 1. We define the threshold density, 9 v (x) by 


6b(x) = 6(5 (x) - ua) 


1 if 5(x) > ua 
0 else. 


( 1 ) 


Note the qualitative difference between 5 which is a weighted density held, and 9 U which 
just defines uniform domains, all having equal weight, and (9 u (x)) = Qi. 

Let us now consider how the biasing changes the distribution in relation to the 
classification we have given in terms of P(0). In what follows we show, for v > v 0 > 0, ( v Q 
given below) 

PM > P( 0) (2) 

where P v (k) and P(k ) are the PS of the biased and underlying held respectively, i.e. 
the Fourier transform of ^(r) and of the normalised underlying correlation function 
£ c (r) = £(r)/<7 2 respectively. This result is independent of £ c (r). The correlation function 
£„(r) of the biased held is given (Kaiser 1984) by the expression 


Qi(v) 2 (Ur) + l) 


1 

2TT^l -£?(r) 




dSdS' 


x exp 


( (6 2 + n-2Ur)55' \ 
{ 2(1 -e c (r)) ) 


(3) 


where the integrand on the right hand side is the two point joint probability density for 
the Gaussian held Using the expression for Qi(v) given above, one can recast this after a 
simple change of variables into the form 


U r ) 


L°° dxe-* 2 ' 2 /; dye~ y2 / 2 
[/“ dxe~* 2 / 2 } 2 


(4) 


where /j = (y — ^ c x )/— £ 2 . In this form it is evident that £ c (r) = 0 £ v (r) = 0 and that 

signer)] = signer)]. 

Taylor expanding this expression about = 0, we find 1 


Zv(r) = h (v)£ c (r) + b 2 (u)^ 2 (r) + ... 


( 5 ) 


1 For the expansion to all orders see Jensen & Szalay (1986). 
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with 

( 6 ) 
(7) 

The first term gives the linear relation obtained by Kaiser (1984), as b\{v) « v 2 for v S> 1, 
valid in the regime |£ c | <C 1 and |^| <C 1. It is easy to check that b 2 {v) is positive definite 
for v > 0 (and b 2 (is) ~ z/ 4 /2 for z/ S> 1), so that to this order in £ c (r) one has the bound 

6 ,( r ) > bi(v)£ c (r) for &,(r) 7 ^ 0 and v > 0 . (8) 

If |&(r)| C 1 at all r this bound suffices to give the desired result (2) for all values of v 
such that b\{v) > 1 . As b\{v) is a monotonically increasing function of u (with 61 ( 0 ) = 2 / 7 r) 
this is equivalent to the requirement v > u a with u Q such that b\{y 0 ) = 1 (i.e. v 0 ~ 0.303). 
To show that there is a value of v above which Eq. (2) is indeed satisfied for all permitted 
values of £ c , it suffices to find the threshold value U\ > v Q such that for all u > V\ one has 

(v ^ 

sign J, t ~bi(u) = sign[Q Mr. (9) 

n^c( r ) 

In fact this is a sufficient condition to have the exact curve £j,(£ c ), given by Eq. (4), all 
above the line £„ = for all ^ c . We have found numerically, using Eq. (4), that this 
condition is satisfied for u\ ~ 0.38. 

Note that, if the condition ( 8 ) holds, this means simply that, relative to the 
asymptotic (|£ c | -C 1 and |^| -C 1 ) linearly biased regime in which £„ « bi(v)£ c , the 
anti-correlated regions are less amplihed (|^| < &i( z/ )l£cl) than the positively correlated 
regions (|^| > 6 i(z/)|^ c |). Thus the integral over the biased correlation function is always 
positive, and the bound (2) thus holds. Further it is easy to see that P„( 0) is finite if P( 0) 
is: is bounded for any value of u, and, has the same convergence properties as £ c at large 

distances. This implies that, if the integral of £ c over all space converges, then also that of 
does. 

In terms of the classification of distributions by P( 0) we thus draw the following 
conclusion: Both the critical type system (with P(0) = 00) and Poisson type system 
(with P(0) = constant > 0) remain in the same class; the super-homogeneous distribution 
(with P(0) = 0) however becomes Poissonian (P I/ (0) = constant > 0). The essential 
reason for these changes is simple: as discussed above the behaviour of the PS is the 
same as that of the mass variance at asymptotically large scales. The biasing process 


h(v) = e" 


,2/2 /r dxxe x /2 

[/“ dxe -* 2 /' 2 ] 2 

rlr.tr . 2 — 1 \p~ P/ 2 


1 _„ 2 /2 f™ dx(x 2 - l)e x ' 

2 VC [f°°dxe -* 2 / 2 ] 2 
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is stochastic in nature, and introduces a variance in the number of objects which is 
proportional to the volume. This new variance will dominate asymptotically over that 
of the original distribution only if the latter is super-homogenous (i.e. its asymptotic 
normalised variance is sub-poissonian, decaying faster than Poisson). Consider for example 
the case of a perfect lattice, which is a super-homogeneous distribution (P(0) = 0) in which 
the normalised variance cr 2 (R) = ((AM(i?)) 2 )/(M(i?)) 2 in a sphere of radius R decays 
asymptotically as 1/R 4 . The distribution obtained by keeping (or rejecting) each point 
with probability p (or 1 — p) is described by a simple binomial distribution, with a variance 
a 2 (R) oc p( 1 — p)/N oc 1 /R 3 ( N being the mean number of points inside a sphere). Biasing 
is not such a purely random sampling, but the effect of stochasticity as a source of Poisson 
variance at large scales is similar. Translated in terms of the PS it gives the result we have 
derived. 

Now let us turn to the implications of this result for cosmological models. Since such 
models have P(0) = 0 in the full matter spectrum, it is evident that we cannot have the 
behaviour P v {k) oc b\(v)P(k) for small k which one might naively infer from the fact that 
£„(r) bi(v)£ c (r) for large separations. Inevitably a non-linear distortion of the biased 
PS at small k relative to the underlying one is induced. How important can the effect be 
qualitatively for a realistic cosmological model? To answer this question we consider the 
simple model PS P(k) = Ake~ k ^ kc . The differences with a cold dark matter (CDM) model 
- which has the same linear Harrison-Zeldovich form at small k but a different (power-law) 
functional form for large k - are not fundamental here, and this PS allows us to calculate 
the correlation function £ c analytically (see GJSL02). This greatly simplifies our numerical 
calculation of the biased PS P u (k), which we do by direct integration of £ v (r) calculated 
using the approximation 

Ur)= \rv||y exp ( 1 ' 2 TTi :)^ 1 (1+0(1 '" I)) (10) 

which is very accurate over most of the range of the integration 2 . In Figure 1 we show P u {k ) 
for various values of the threshold v — 1, 2, 3. We see that the shape of the PS at small k is 
completely changed with respect to the underlying PS. Indeed the main feature of the latter 
in this range - the display of a clear maximum and “turn-over” - is completely modified. 

2 This approximation is obtained by expanding the full expression for £„( r ) given in Eq. (4) in and 
further assuming only that V\J (1 — £ c )/(I + £ c ) ^>1. It is a much better approximation than that of Politzer 
and Wise both at small and larger values of £ c . In particular it gives an asymptotic behaviour £„ w ( v 2 + l)£ c 
for u 2 £, c “C 1 which is a much better approximation to the exact behaviour at typically relevant values of v 
(6i(l)«2.4) 
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Qualitatively it is not difficult to understand why this is so. The only characteristic scale 
in the PS (and also in the correlation function) is given by the turn-over (specified in our 
case by k — k c ). On the other hand, the value of P u ( 0) is just the integral over all space 
of £„ which is proportional to the overall normalisation A and (since it is strictly positive) 
must be given on dimensional grounds by Ak c times some function which depends on v. 
For v ~ 1 this function is of order one, so that P v ( 0) ~ max[P(k)]. 

This last point is better illustrated by considering the integral Ja(r, v) = 47r / 0 r x 2 £, u (x)dx 
which converges to P u ( 0) = lirn^oo J 3 (r,u). In Figure 2 the value obtained for it by 
numerical integration of the exact expression given by Eq. (4) for £ v (r) is shown for 
u — 1,2, 3. We also show the same integral for £ c which converges to P( 0) = 0. While 
the latter decreases at large r, converging very slowly to zero (as 1 jr since £ c (r) oc — 1 /r 4 
at large scales), the former all converge towards a constant non-zero value We see that 
the integral picks up its dominant contribution from scales around (and above for v = 1 ) 
r ~ 10 Mpc (see caption for explanation of the normalisations, which are irrelevant for the 
present considerations). From the inset in the figure, which shows both £ c (r) and £j,(r), 
we see that this is the scale below which the correlation function is non-linearly amplified. 
Moreover it is shown that the smaller scales at which £„(r) is most distorted relative to 
£ c (r) do not contribute significantly to J 3 (because of the r 2 factor). This fact also explains 
the accuracy of the PS obtained using the approximation Eq. (10) for £„(r), which can 
be seen by comparing the asymptotic values of the integrals in Figure 2 with P l/ ( 0) in 
Figure 1. Note that for v — 1 the distortion away from linear is relatively weak in the 
part of the correlation function which dominates the integral in J 3 (r, u ), and that there is 
even a non-negligible contribution from the larger scales at which the correlation function 
amplification is extremely close to linear. 

Let us now draw our conclusions. We have calculated the effect of the distortion 
at small k of a biased PS relative to an underlying PS, which we have shown to be an 
inevitable effect of biasing on cosmological models. In particular we have used a simplified 
model PS and the simplest (but reference) threshold biasing scheme of Kaiser. The latter is 
not a realistic model of biasing for various reasons, most notably because of the extremely 
strong non-linear distortion of the two-point correlation function at small scales whilst the 
observed ratio of correlation between galaxies and clusters is approximately linear. Our 
results however should be qualitatively correct for any biasing scheme and standard CDM 
type model. The only thing which we expect to depend on the biasing scheme is the 
functional form of the PS to the right of the turn-over (for k > k c ), and correspondingly the 
shape of the correlation function at small scales, which will only make a minor numerical 
change to our calculation. The essential feature of biasing which brings about the effect we 
have discussed - a non-linear amplification of the correlation function - will be common to 



any biasing model. In fact the converse of what we have argued is that any biasing scheme 
which is stochastic (i.e. any procedure for selecting sites for objects which is probabilistic) 
must give such a non-linear distortion of the correlation function for cosmological models: 
if the correlation function were amplified linearly at all scales, the asymptotic behaviour of 
the variance will be sub-Poissonian instead of Poissonian. 

We finally draw two conclusions with respect to the comparison of cosmological models 
with observational data. In order to make the link to observations of galaxies or clusters, 
current theories make use normally of biasing in one form or another. Invariably however 
the measured PS from observations is fit to a model spectrum which is just a rescaled dark 
matter PS (for a recent example see Lahav et al. 2002). Our first conclusion is that such a 
behaviour cannot be obtained by biasing as the PS is not linearly amplified neither at small 
nor at large wavenumbers. Our second conclusion is that it will be very useful to measure 
not just the PS at small k, but also the real space correlation function at large distances. 
Invariably only the first is considered in analysis of observations of galaxy distributions at 
large scales, because, it is argued, it is the natural probe given that cosmological models of 
structure formation are formulated in k space. We have seen however that biasing leads to a 
distortion of the underlying dark matter PS, while (insert panel of Figure 2) the correlation 
function at sufficiently large distance remains undistorted. In particular for standard models 
with a H-Z PS at small k, the cleanest way to detect this behaviour in biased objects is 
by looking at the correlation function at large scales, which should maintain the behaviour 
£ u (r) ~ —r” 4 associated with the small k behaviour of the underlying dark matter PS. To 
address the viability of measuring this behaviour in current and forthcoming surveys of 
large scale structure requires in particular the detailed treatment of both the passage to 
the discrete distribution (we have treated here always continuous density fields), and the 
question of the variance of estimators of £(r). 
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Fig. 1.— The PS P v {k ) derived from the biased correlation function £ v (r) for values of 
the threshold v — 1,2,3 is shown. The underlying correlation function which gives £ c (r) is 
that derived from the P(k) = Ake~ k ^ kc and the approximation given in Eq. (10) for £ v is 
used. The clear distortion of the PS at small k is seen, the “turn-over” in the underlying 
PS essentially disappears already for u — 1. The constants k c and A are fixed by £ c (0) = 1 
and the requirement that £ c (r) = 0 at r = 38 Mpc. The latter is taken as in a typical CDM 
model (see e.g. Padmanabhan 1993). We could alternatively fix the wavenumber at the 
maximum of the PS. 


Fig. 2.— The integral J 3 (r) for the same underlying correlation function as in Figure 1 
and for the same range of values of u, calculated with the exact expression Eq. (4) for 
£„(r). Also shown is the analagous integral of £ c . While the latter converges slowly to its 
asymptotic value of P( 0) = 0, the other integrals converge to constant non-zero P v { 0). They 
are dominated by the range r ~ 10 Mpc where, as can be seen from the inset which shows 
both £ c and each of the £„(r), the correlation functions £„(r) are amplified non-linearly. The 
contribution from the extremely amplified region at small r is small (because of the r 2 factor 
in the integral), which also makes the approximation in calculating the PS with Eq. (10) 
very accurate, as can be checked by comparing the asymptotic values with those of P u ( 0) in 
Figure 1. 



